Post-processing with linear optics for improving the quality of single-photon sources 
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■^j- i Triggered single-photon sources produce the vacuum state with non-negligible probability, but 

produce a much smaller multiphoton component. It is therefore reasonable to approximate the 
output of these photon sources as a mixture of the vacuum and single-photon states. We show that 
' it is impossible to increase the probability for a single photon using linear optics and photodetection 

on fewer than four modes. This impossibility is due to the incoherence of the inputs; if the inputs 
were pure-state superpositions, it would be possible to obtain a perfect single-photon output. In 
the more general case, a chain of beam splitters can be used to increase the probability for a single 
photon, but at the expense of adding an additional multiphoton component. This improvement is 
robust against detector inefficiencies, but is degraded by dark counts or multiphoton components in 
the input. 
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I. INTRODUCTION 



>> ! PACS numbers: 03.67.-a, 42.50.Dv 
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One of the most promising methods for quantum information processing is linear optics and photodetection. Linear 
optics and photodetection may be used for provably secure quantum communication p|, as well as quantum compu- 
tation 0. An important requirement for these schemes is the ability to produce a single photon on demand 
Q_i! yet generating high-fidelity single-photon states is challenging. The traditional method for generating single photons 
involves photodetection on one output mode from a non-degenerate parametric down-conversion process to post-select 
a single photon in the correlated mode 0, |B| . This method has the drawback that the time of the photon emission is 
not controlled. More recently, triggered photon sources have been developed, including molecules [g, quantum wells 
0, colour centers Q, ions and quantum dots Hfl. These sources have a significant vacuum contribution, but the 
multiphoton contribution may be made very small 1 1 1| . 

For the majority of this study, we approximate these sources by taking the multiphoton probability to be zero. 
That is, we consider an idealised single-mode single-photon source, which may be represented by the density operator 



p p = (l-p)|0><0|+p|l)<l|. (1) 

Here p is the probability for a single photon, and is also called the efficiency. Increasing the efficiency is important 
because many quantum optics experiments, especially those concerned with linear optical quantum computation, 
require high efficiency sources. We also derive results for states with a coherent superposition of the vacuum and a 
single photon; however, this form of state is not produced by single-photon sources. 

Much effort is directed towards improving sources, but here we pose the question as to whether it is possible to 
perform post-processing to obtain higher efficiency. Ideally this post-processing should also maintain a zero multi- 
photon contribution (though a very small multiphoton contribution would also be acceptable). A promising method 
of post-processing is linear optics and photodetection. As mentioned above, linear optics and photodetection can be 
used to perform quantum computation and optical controlled-NOT gates have recently been demonstrated [l2j : 
however, there are also no-go theorems for linear optics [l3|. 

In a recent publication [l4j we investigated the possibility of improving the single-photon efficiency p with linear 
optical elements and photodetection. Here we present a number of new results that clarify the limitations inherent 
in this method, as well as reviewing the results presented in [T^. In particular, we analyse the effects of various 
experimental limitations, present a scheme that gives perfect results for inputs with a coherent superposition of zero 
and one photon, and show connections between some difficult unsolved problems. 



2 



We begin by showing that it is not possible to obtain an improvement for the simple case of a beam splitter 
in Sec. \n\ In Sec. 11111 we show that, if we allow a coherent superposition of zero and one photon, rather than an 
incoherent mixture, it is possible to obtain a perfect single-photon state. We then proceed to the case of a multimode 
interferometer with incoherent inputs in Sec. IIVI In Sees. and IVII we give an expanded discussion of the limit 
on the improvement that it is possible to obtain, and the method to obtain an improvement. Sec. IVIII gives a 
detailed discussion of the impact of various experimental problems on this method. We give further discussion of 
no-go theorems for post-processing in Sec. IVIIII Then, in Sec. IIXI we show that there are deep connections between 
the unsolved problems for post-processing. Lastly we give a discussion of how the theory is changed by allowing 
multiphoton contributions in the inputs in Sec. an d we conclude in Sec. IXII 



II. BEAM SPLITTER 



To begin, we consider the simplest case of two copies of the quantum state combined on a single beam splitter. 
The initial state can be written in the form 

P?J = (1 -Pi)(l -p 2 )|00)(00| +pi(1-P2)|10><10| + P2 {1 -pi)|01)(01| + p lP2 |ll)(ll|. (2) 

Each photonic mode operator gets transformed by the beam splitter in the following way: 

a\ i-> A n a\ + A 2 i<4, a\ h-> A 12 a\ + A 22 a2,, (3) 

where A is a 2 x 2 unitary matrix |15|. We have two options, to project onto the vacuum or onto the single-photon 
Fock state (projecting onto two photons results in vacuum output in the unmeasured mode). After vacuum projection 
in mode 2, we end up with an unnormalised state of the form 

*P0 ~ ln\/nl -L ( Pl i a 12 



PXt oc |0)(0| + j^\A n \ 2 + ~^\A 12 \ 2 |1)(1| + 2-photon term. (4) 

V 1 - pi 1 - vi J 

Thus the ratio between the probabilities for one and zero photons is just a weighted average of Pi/(1 — Pi) and 
P2/(l — P2), and cannot exceed either of these. That is, it is not possible to improve the ratio between the probabilities 
for obtaining one and zero photons. This automatically implies that it is not possible to improve the absolute 
probability of obtaining one photon. 

In the case of a single-photon detection, the resulting unnormalised state is 

J2) _ fl-pi,. |2 l-p 2|A ,2^ ,r,\/nl , I A |2| 



iCt « — ^|A 22 | 2 + -\A 21 \ 2 |0)(0| + | P erAHl><l|, (5) 

V Pl P2 J 

where perA is the permanent |16| of the beam splitter matrix A. Since the absolute value of the permanent of a 
unitary matrix is bounded from above by unity, and the term in brackets is a weighted sum of terms (1 —Pi)/pi, we do 
not find any improvement for this case either. Hence there is no improvement in the probability for a single photon 
if zero, one or two photons are detected. These results demonstrate that for mixed-state inputs it is impossible to 
obtain an improvement in the single-photon probability using a beam splitter. 

III. PURE-STATE INPUTS 

It is possible to obtain an improvement using a beam splitter if the inputs are in pure-state superpositions of zero 
and one photon, instead of incoherent mixtures. Consider two input modes that are each in the state a|0) + 
The initial state may be written as 

Win = [« 2 + + 4) + P 2 44} |00). (6) 

Applying the beam splitter transformation J2J| gives 

Wir-L = {a 2 + [(An + A 12 )a\ + (A 21 + A 22 )a 2 ] + p 2 [AnAi 2 (a t 1 ) 2 + A 2 iA 22 (a 2 ) 2 + (perA)<44] } |00). (7) 
Conditioning on detection of zero photons in mode 2 gives the output state 

IV^t (x a 2 |0) + a/3(A n + A 12 )|l) + \/2> 2 A u A 12 |2). (8) 
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It is easily seen that this output state may have a higher probability for a single photon. For example, if the initial 
state is close to the vacuum state (i.e. a > then an improvement by a factor of two may be obtained by using 
An = Aia = 1/V2- 

In fact, it is possible to further process this output state to obtain a perfect single-photon state. If we combine this 
state with mode 3, which is also assumed to be prepared in state a|0) + /3|1), then the total state may be represented 
as 



v(3) 



oc 



a 



■ ap(A n + A 12 )a{ + /3 2 AnAi 2 (a t 1 ) 2 ] (a + /3aj) |00). 



(9) 



Applying the beam splitter transformation J3} (except using a prime to distinguish this beam splitter from the previous 
one) gives 



(3) 



OC 



a - 



/3(Ai 3 al+A^a|)l|0Q). (10) 



+ a/3(An + Ai 2 )(A' 11 a t 1 + A 31 a|) + /3 2 AnAi 2 (A' 11 a t 1 + A^aJ,) 2 
Conditioning on detection of two photons gives 

IV^t « V2f3 2 A' 31 {a [(An + A 12 )A' 33 + A n A 12 A 31 ] |0) + /3A u Ai 2 (2A' U A 33 + A' 31 A' 13 ) |1)} 
To make this equation more clear, we use 



A = 



e' v cos ( 
sin 



smtf 
^ cos ( 



A' = 



e^' cos ( 
sin0' 



sin 0' 



cosf 



The condition that the output state is a pure single-photon state then becomes 

(cos - e" 10 sin 0)e"^' cos 9' - sin 9 cos 9 sin 9' = 0. 
This equation may be satisfied by taking 9' and <f>' to be 



arctan 



sin e/ cos I 



arg(cos 9 — e ^ sin 



(11) 



(12) 



(13) 



(14) 



That is, regardless of the characteristics of the initial beam splitter (provided sin 9 and cos 9 are nonzero), it is possible 
to obtain a perfect single-photon output. 

Another issue is the probability for obtaining the desired pattern of detection results. Using the unnormalised 

(3) 

expression for |V7out above, this probability is given by 



sin 



2 Q/\2 



P = 2| / 3| 6 sin 2 0cos 2 0sin 2 0'(2cos 2 

^gisin 2 29{l - cos</> sin 29) (^sin 2 20- 1 + cos <j> sin 26 
(isin 2 20 + 1 -cos </>sin 20) 3 



(15) 



In the second line we have used the expression (|14H for 0'. This probability is plotted in Fig. U for the range to 7r in 
and 4> {P is periodic with period 7r in these variables). There are four maxima in this range, for (0, </>) = (n/i, tt), 
(7r/4,acos(13/14)), (3tt/4,0) and (3tt/4, acos(- 13/14)). The exact values of </> of acos(13/14) and acos(-13/14) are 
not obvious from the plot but are straightforward to obtain analytically. The two maxima (7r/4,7r) and (37r/4, 0) 
correspond to the same beam splitter, so there are only three maxima that correspond to distinct beam splitters. 
Each of these maxima is exactly the same height, 



= 16|/3| 6 /81. 



(16) 



One factor that distinguishes the maxima is the sensitivity to the parameters. Clearly the maxima at 
(7r/4, acos(13/14)) and (37r/4, acos(— 13/14)) are far more sensitive to the values of and 0, and it is therefore 
better to use the beam splitter corresponding to (7r/4, tt) and (37r/4, 0). For the second beam splitter, the appropriate 
parameters are 0' = acos(l/3) and (/»' = 0. That is, the best result is obtained by using a 50/50 beam splitter followed 
by a beam splitter with a reflectivity of 1/9. 

Thus we see that, if the inputs to an interferometer are in pure superposition states, it is possible to obtain a perfect 
single-photon output for three modes. In contrast, if the inputs to the interferometer are incoherent superpositions 
of Fock states, it is impossible to obtain an improvement in the single-photon probability for three modes [lj|- These 
results imply that it is the incoherence in the inputs that prevents an improvement in the single photon probability. 
It would be interesting to determine the degree of decoherence that is sufficient to prevent an improvement in the 
single-photon efficiency. However, this is a difficult problem, which we leave to further investigation. 
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FIG. 1: The probability of obtaining the desired detection results as a function of the beam splitter parameters 9 and (f> for 
the first beam splitter. 
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FIG. 2: Schematic setup of the network. We assume -/V incoming modes prepared in the state $17\ with different pi. The 
photon number is measured in output modes 2 to N, and we wish to improve the probability for a single photon in mode 1. 

IV. MULTIMODE INCOHERENT INPUTS 

Although it is possible to obtain perfect single-photon states from pure superposition states, this method can not 
be applied to current experiments, as single-photon sources do not produce pure superposition states. Real single- 
photon sources produce an incoherent combination of Fock states; therefore we consider input states of this form for 
the remainder of this paper. In the multimode case we start with a supply of N mixed states of the form l[TJl. For 
additional generality we allow the different inputs to have different probabilities for a single photon, pi, and we denote 
the maximum of these probabilities by p max . The initial input state may be described by 

N 

A ( „ w) = (8)[(i-^)|o)(o|+ P4 |i>(i|] 

i=l 

=E p «(n( a i) si i o )<°in^) si )' 

s \ i i / 

where P s = YliP? (1 — Pi) 1_Si > an d the vector s — (s\, • • • , sn) t , {si = 0, 1), gives the photon numbers in the inputs. 
The quantity P s is the probability of obtaining this combination of input photon numbers. 

This input is then passed through a passive interferometer which consists of beam splitters, mirrors, and phase 
shifters. Each of these elements preserves total photon number from input to output under ideal conditions. No 
energy is required to operate these optical elements, hence the term passive (also known as linear optical elements). 
More generally, polarisation transforming elements can be included, but here we are concerned only with a scalar 
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field treatment; in fact polarisation effects could be included by doubling the number of channels and treating the 
two polarisations in a mode as two separate channels. 

Classically the field amplitude of channel i would be represented by the complex number a,-. The set of all field 
amplitudes for the TV-channel interferometer is given by the vector a = (aj., • • • , aAr) T . The passive interferometer 
transforms the input amplitudes to the output amplitudes via the matrix transformation a i— > A'a with A £ U(iV), 
where U(iV) is the set of all N x N unitary matrices. Quantisation of the field is obtained by the replacement of a 
by the vector annihilation operator a, and the interferometer transforms the operators according to a} i— > A T d^ 
This transformation of the operators yields 



r trans 




(18) 



In the completely general case, we could perform photodetections on N — Ni of the modes, and use the remaining N% 
modes as single-photon sources if the desired combination of detection results is obtained. However, this generality 
is not needed here because we are concerned with the maximum improvement in the photon statistics in a single 
mode. In order to fix notation, we denote as mode 1 that mode for which we want to improve the statistics, and label 
the other modes where photodetections have not been performed as modes 2 to N\. The reduced density matrix in 
mode 1 is then identical to what would be obtained if photodetections were performed on modes 2 to Ni (as well as 
Ni + 1 to N), and the results of these photodetections discarded. Therefore, the probability for a single photon will 
be a weighted average of the single-photon probabilities for the different combinations of detections in modes 2 to 
Ni . Hence the maximum single-photon probability in mode 1 will be obtained for some combination of detections in 
modes 2 to iVi . For this reason we consider the state in mode 1 conditioned on photodetections in the other N — 1 
modes. As our aim is to determine the best results possible using linear optics and photodetection, we also assume 
that the photodetectors perform perfect photon counting measurements (imperfect detection is discussed later). 

Before determining the conditional output state, we introduce some additional notation. The total number of 
photons detected is D, and the maximum possible number of photons input to the interferometer is M . As some of 
the pi may be equal to zero, M may be less than N; M is equal to the number of nonzero values of pi. For j > 1, 
rij is the number of photons detected in mode j , and n\ is the photon number in mode 1 (the output mode). We use 



the notation E n = n, (so £„ = D + n±) and £ s = J^. s,-. In addition, we define the set <I> S 
Y s be the set that consists of all vectors comprised of the elements of (f> s . 
The conditional state in mode 1 after photodetection in modes 2 to N is 



{i\si = 1}, and let 



Each coefficient c ni is given by 



JV 

Pont = X! C »iK)("ll 

711=0 



c ni = K(n\p[^J as \n), 



(19) 



(20) 



where \n) is a tensor product of number states in each of the output modes and the normalisation constant K is equal 
to 



K 



Evaluating c ni gives 



Yl w i-i"' 

.ni=0 



K' 



nv. 



^ ] Pa \Sa,n\ j 



(21) 



(22) 



where K' = Kj n j*> ano ^ 



S a 



E (Ai,^-Ai,0-(A 



<t6Y, 



-n N + l 



This quantity may alternatively be expressed using permanents as 

S s ,n = per(A[n, s]). 



•••A 



N,<rs, 



(23) 



(24) 
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Here the notation A[n, s] is used to indicate that the i'th column of A is repeated s, times, and the j'th row is 
repeated rij times. 

Two figures of merit for an arbitrary single mode field cft|i)(i| are 

R=^, G = ^i. (25) 

qo qi /qo 

We use the subscript "out" to indicate the output field, and "in" to indicate the input field. For the output field we 
simply have qi = Ci. For the input field, we have G m — 0, as the two-photon component is assumed be negligible. For 
simplicity we define i?j n to be the maximum input ratio p m ax/(l — Pmax)- 

The figure of merit G characterises the two-photon contribution, and is equal to 1/2 for Poisson photon statistics. 
If the multiphoton component in the output is zero, then comparing i?j n and i? ut immediately tells us if there is an 
improvement in the probability for a single photon. Even if the multiphoton component is nonzero, using i? ut has 
the advantages: 

1. The common constant K' cancels, so it is possible to evaluate i? ut analytically. 

2. If Rout < R'm, then it is clear that c\ < p max - Thus we can determine those cases where there is no improvement. 

3. For p max <C 1, Co ~ 1 and i?i n ~ p max - Therefore the improvement in R is approximately the same as the 
improvement in the single-photon probability over p max - 

An alternative measure of the multiphoton contributions is given by how sub-Poissonian the field is. That is, we 
may define the measure 

n=<^. (26) 

(n) 

For a sub-Poissonian field, II < 1. States of the form are sub-Poissonian, with II = 1 — p. We take IIj n to be the 
minimum value in the inputs, 1 — p max , and n out is simply the value for the output mode. If an output c\ greater 
than Pmax is obtained, while maintaining a multiphoton contribution that is zero or very small, then it is clear that 
n out will be smaller than IIi n . On the other hand, if the output has multiphoton contributions similar to those for a 
Poisson distribution, n out will be closer to 1. 



V. LIMIT ON IMPROVEMENT 



Ideally we wish to obtain an improvement in the figure of merit R, while maintaining a value of G that is zero, or 
at least small with respect to 1/2 (the value for a Poisson distribution). This is a difficult task, so for simplicity we 
begin by focusing on improving R. As was shown in Ref. |14| . there is an upper limit on how far R can be increased. 
Here we show this result in more detail. 

First we consider the expression for Cq: 

C = K' P s\ S s,n,°\ 2 - ( 27 ) 

s;S„=D 

Here n® = and n® is the combination of detection results for j > 1. A simple way of re-expressing this summation 
is 

c = N _ D ^ Ps k \S s >*,n°\ , (28) 

s;S s =D+l k;s k =l 

where sf = Si except for s\ = 0. That is, we consider combinations of input photons s with one too many photons, 
then remove one of these photons to obtain the correct number of input photons. This expression for the sum gives 
each term N — D times, so it is necessary to divide by N~ D to obtain the correct result. Specifically, each alternative 
s k may be obtained from an s which is identical, except one of the zeros of s k is replaced with a one. As each s k has 
N — D zeros, there are N — D possible alternative s that give the same s k . 

If some of the inputs to the interferometer are simply vacuum states (i.e. some of the pi are zero), it is possible to 
express the summation in a more efficient way. First note that those terms in the sum in Eq. (|27(l with P s = do not 
contribute to the sum; therefore we may restrict to terms with P s ^ 0. 

c Q =K' J2 P *\ S s, n °\ 2 - (29) 

s:P„=iO 
T. a =D 
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We may re-express this equation as 



K' 



Co 



M-D 



E E p s» I* 



(30) 



s;P fl ^0 fe;s fc = l 



where = S a k n o. Here the dividing factor required is only M — D. 

Recall that the maximum total number of photons is M, so there are N — M inputs with pi = 0. Each alternative 
s k still has N — D zeros, but some of these zeros will correspond to inputs with pi — 0. Because we are restricting 
to terms with P s ^ 0, for all i with pi — 0, Sj = and therefore s\ must be equal to zero. Thus all N — M of the 
inputs with pi = must correspond to zeros of s k , and so there will only be M — D zeros of s k that correspond to 
nonzero Pi . As in the previous case, each s k may be obtained from an s which is identical, except one of the zeros of 
s k is replaced with a one. However, because we have restricted to terms with P s ^ 0, the zero that is replaced with 
a one must be for an i with nonzero pi. Hence there are only M — D alternative s that lead to the same s k , and the 
redundancy in this case is only M — D. That is why a dividing factor of M — D is required in this case. 

Simplifying Eq. I)3()[l . wc obtain 



K' 



Co 



E Ps E — i^-i 2 



M-D ^ <—> Pk 

T. B =D+l 



(31) 



Because the sum is limited to terms where P s ^ 0, pk is nonzero, and therefore the ratio (l—pk)/pk does not diverge. 
Using pk < Pmax, we obtain the inequality 



c > 



K'/R h 
M-D 



E p- E i^i 2 - 



(32) 



s;S B =D+l k;s k = l 



We now allow terms with P„ = 0, because they do not contribute to the sum. 
It is also possible to obtain an inequality for c± . The probability c\ is given by 



Cl 



K' J2 P ° \ S « 



(33) 



s;S B = D + l 



In this case the notation n 1 means n\ = 1, and nj is the combination of detection results for j > 1. We may express 

Sen. 1 3-S 
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s.** 1 ~~ E 



Ai, CTl (A 2i(T2 • ■ ■ A 2 , CT „ 2+1 ) ■•• f 

E A ^ E ( a 2, ct1 ---a 2 ,.„ 2 )- 

E ^lkSs k ,n° 
fe;s fc = l 

E A i^fe- 



A 



■ A 



A 



AT.cts 



fc -"JV + 1 



A 



JV,<T£ 



fc;Sfc — 1 

Therefore, we may re-express the equation for c\ as 



d=K' p° 



E A-ikSk 



k:s k =l 



This then gives the inequality 



Cl 



E p * E i^i 5 



Combining Eqs. 1)32(1 and (|36[1 . we can see that 



;£.=D+1 fe:si. = l 



Cl 



flout = — < Rin(M-D). 
CO 



(34) 



(35) 



(36) 



(37) 
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This yields an upper limit on the ratio between the probabilities for one and zero photons. This result allows one to 
draw three main conclusions: 

1. As M < N and D > 0, the improvement in R can never be greater than TV. There is no known scheme that 
saturates this upper bound, but there is a scheme known that achieves an improvement of approximately N/4 [T^ |. 

2. If the number of photons detected is one less than the maximum input number, then M — D = 1, and there 
can not be an improvement 19J. This case is important because it is the most straightforward way of eliminating 
the possibility of two or more photons in the output mode. We have not proven that it is impossible to obtain an 
improvement while eliminating the multiphoton component, but if such a scheme is possible it can not eliminate the 
multiphoton component by detecting one less than the maximum input number of photons. 

3. It is impossible to obtain a single photon with unit probability if p max < 1. If Ci = 1 were obtained, then i? ou t 
would be infinite; from Eq. (|37[) . this is clearly not possible unless i?; n is infinite (which would correspond to p max = 1). 



VI. METHOD FOR IMPROVEMENT 

We have previously shown that it is possible to obtain an improvement in the probability for a single photon |14| . 
Here we review this method, giving more motivation for this scheme and propose a simple realisation using a line 
of beam splitters. In order to obtain a value for the ratio i? ut that is close to the upper limit, we require the two 
inequalities (|32[1 and (|36fl to be as close to equality as possible. We may achieve equality in the first case (|32[l by 
taking all nonzero pi equal to p max . 

To obtain equality in Eq. I|36|L we would require s/. = 1 whenever Ai/. is nonzero, and Sk to be proportional to A^ fe 
for those values of k where Sk = 1. These conditions would need to be satisfied for all s that give nonzero P s . The 
first condition is a problem, because it can not be satisfied unless D = M — 1. To see this, note that this condition is 
equivalent to requiring that Sk = implies Aifc = 0, for all s that give P s > 0. If P s > 0, then Sk = for all k such 
that pt = 0. Therefore, pk = implies A\h = 0. Now if D + 1 < M, then it must be the case that Sk — for some k 
such that pk > 0. In addition, for arbitrary k, there will be an s with P s > such that Sk — and pk > 0. Therefore, 
for the first condition to be satisfied, it would be necessary for Aik to be equal to zero for all k. This is clearly not 
possible, because A is a unitary matrix. 

On the other hand, if D + 1 = M, then the only s giving P s > is that with Sk — 1 for pk > 0, and Sk = for 
Pk = 0. Therefore, it is possible for the first condition to be satisfied, by choosing a A with Axk = for all k such that 
Pk = 0. However, the case with D = M — 1 is unimportant, because it is not possible to obtain an improvement in 
R. Hence we see that, in any case where it is possible to obtain an improvement, it is not possible to obtain equality 
inEq. 

On the other hand, we can determine a scheme that gives Sk oc A* fc . This condition can be satisfied using the 
interferometer with matrix elements 

An = -e, A 2 i = v 7 ! - e 2 , 
Ah = v/(l-e 2 )/(/V-l), k 2l = e/VN~T, (38) 

for i > 1 (the values of A^ for i > 2 do not enter into the analysis). Here e is a small number, and we ignore terms 
of order e or higher. Now let pi = p max , and consider the measurement record where zero photons are detected in 
modes 3 to N, so the number of photons detected in mode 2 is D. 

This scheme is the same as in Ref. ^4|j except input modes 1 and 2 have been swapped. Expressing the scheme in 
this form allows us to determine a simple realisation using beam splitters. This scheme may be performed using the 
chain of beam splitters shown in Fig. El The first A — 2 beam splitters (those on the right) result in an output beam 
with equal contributions from N — 1 of the inputs. This equal combination is achieved by decreasing the reflectivities 
from 1/2 for the first beam splitter to 1/(N — 1) for beam splitter N — 2. The last beam splitter has the low reflectivity 
e 2 . With appropriate phase shifts, these beam splitters give the overall interferometer described by A in Eq. (J3SJ. 

To determine c ni , note first that A21 ^> A 2 i for i > 1, so we may ignore those terms in the sum for S s . n where A21 
does not appear. Each term has magnitude A™2 A 2 iA 2 ^~ 1 [20|. and there are D(D + n\ — 1)! such terms. Therefore, 
provided s\ = 1, 

S., n w D{D + m - lJIAgAaiAg" 1 . (39) 

If si = 0, then S s . n is of order e. 

Recall that Sk — S s k n o, where s k is equal to s, except for s| = 0, and n\ — 0. (We do not consider k where 
Sk = 0.) The result for k > 1 may be obtained by replacing n\ with and D with S s — 1 in Eq. 139fl . giving 

5 , fc «(S s -l)(E s -2)!A 2 iAf ; ?- 2 . (40) 



9 





FIG. 3: A realisation of the interferometer specified in Eq. I|38|l using beam splitters. The reflectivities of the beam splitters 
are specified in the labels on the beam splitters. This example is for N = 6. 



For k = 1, we simply obtain Si of order e. Similarly, is constant for k > 1, and of order e for k = 1. Thus this 
scheme gives Sk oc A* fc , as claimed above. 

In order to determine c ni , note that there are ( D+ ~ _i) different combinations of inputs such that S s = D + n\ 
and si — 1. Combining this expression with Eq. (|39f) . we have 

c ~ K ' n n +ni(1 _ n sn-p-n, (N - l)\D 2 {D + m - 1)! 2n 2 2g _ 2 
m! x x (JV-D-ni)! ^12^21^22 

= jr f_^L,r y + m-D! (41) 
\N-lJ m!(JV-D-ni)! 



We have combined those factors that do not depend on ri\ into a new constant if", and used A 12 ~ ^/yN — 1. Using 
Eq. gives 

D(N — D) 

Rout ~ Rin _ 7 • (42) 

The maximum improvement in i? is obtained for Z? = [iV/2], where i? ut ~ Rin[N 2 / 4\ / (N — 1). The multiplicative 
factor [A^ 2 /4J / (N — 1) is larger than 1 for all N > 4. Thus we find that, provided there are at least 4 modes, we may 
obtain an improvement in i? ou t- For p max < 1, Cj « _Pmax|_-^ 2 /4j/(-^ — !)■ For large N, the probability of a single 
photon increases approximately as AT/4, but does not achieve the upper bound of N. 

Although we find an improvement in the measure R, the two-photon contribution is not negligible. Using the 
measure G, we find 

r - C2 / Cl ~ (D + 1)(N-D-1) 

Gout " ^ 2D(N-D) • (43) 

For D = fiV/2] , this measure is close to 1/2, so the two-photon component is similar to that for a Poisson distribution. 
By taking D = N — 2, it is possible to obtain an improvement in R of about a factor of two, with a value of G ou t 
about half that for a Poisson distribution. However, this two-photon contribution is still much greater than for good 
single-photon sources [llj. 

The multiphoton contributions are especially important for larger p max . Although the improvement in R is inde- 
pendent of Pmax, the multiphoton component means that improvements in ci are obtained only for values of p mg , x 
below 1/2. That is, this method can only be used to obtain improvements in the probability of a single photon up to 
1/2, but not to make the probability of a single photon arbitrarily close to 1. 

This scheme also performs poorly when evaluated via the measure II. It is more difficult to evaluate this scheme 
using this measure, and we do not have a simple solution for II. However, numerically we find that n out is close to 1 
for D = \N/2] , again indicating that the output state is close to Poissonian. For D = N — 2, n ou t is closer to Hi n 
(= 1 — p m ax), but for no D is a value of n out less than Hj n obtained. 

Nevertheless, this scheme does give an improvement in R. It can be expected that this improvement is close to 
the maximum possible, because this scheme satisfies Sk °c AJ fc . However, note that it satisfies the other condition 



10 



for optimality, i.e. that An 7^ implies s& — 1, fairly poorly. As shown above, this condition can not be satisfied 
completely (unless there is no improvement in R), and there does not appear to be any method of satisfying it better 
than the method we have described above. Extensive numerical searches have failed to find any scheme that gives a 
better improvement in R than the above scheme, strongly indicating that it is optimal for increasing R. 

VII. EXPERIMENTAL LIMITATIONS 

In practice, there will be a number of limitations to using this method for improving the probability of a single 
photon. The main ones are: 

1. Real photodetectors do not give perfect photon counting measurements. Most photodetectors can only distinguish 
between the vacuum state and a state with one or more photons. 

2. Real photodetectors have limited efficiency and dark counts. 

3. The desired combination of detection results will occur with low probability. 

4. Real sources will have a finite multiphoton component. 

For N — 2 of the detectors, point 1 will not be a problem. The reason for this is that we are conditioning on 
detection of zero photons at these detectors. It is only the detector on mode 2 that is required to perform a photon 
counting measurement. Even for this detector, it is not necessary to determine the exact photon number. That is 
because the probability of a single photon will be increased for any number of photons detected from 2 to N — 2. 
Therefore, if the detector can register that the photon number is in this range, rather than the exact photon number, 
it will be sufficient to produce an improvement in the single-photon probability. 

Alternatively, an improvement can also be achieved using detection that simply verifies that there is more than one 
photon. For example, the Visible Light Photon Counter can do this task with high efficiency. Even though the 
possibilities of N — 1 or N photons have not been eliminated, they have lower probability, and will not contribute 
significantly to the photon probabilities. 

Limited efficiency is not a severe problem for small p ma x, because the probability for the vacuum is relatively large. 
Dark counts will not be a problem for the first N — 2 detectors, because we are conditioning on vacuum detection 
at these detectors. Dark counts will merely slightly reduce the probability of obtaining the desired combination of 
detection results. On the other hand, dark counts will be a problem for detector 2, as we are conditioning on detection 
of more than one photon at this detector. 

Point 3 will always be a problem, because the above scheme is only effective for small e. The probability for this 
combination of detection results becomes very small in the limit of small e. If larger values of e are used, then the 
final probability for a single photon becomes smaller. Thus there is a trade-off involved. Point 4 is more complicated, 
and it is not clear how important this problem is without performing direct calculations. 

To estimate the relative importance of each of the above problems, we have calculated the final conditional prob- 
ability for a single photon successively taking each of the above issues into account. In Fig. we have plotted the 
conditional probability for a single photon versus the probability for obtaining the desired combination of detection 
results. These curves are parametrised by e; that is, both probabilities were calculated for a range of values of e. In 
general, as e is decreased, the probability for obtaining the desired detection results decreases, and the final conditional 
probability for a single photon increases. The particular example we have shown is of a 4 mode interferometer where 
Pmax = 0.2 for the inputs. 

To perform these calculations, the density matrix was left unnormalised. The trace of the density matrix at the 
end of the calculation then gives the probability for that combination of detection results. For those cases where 
two, three and four photons were not distinguished, the density matrices for these three cases were simply added. To 
take account of finite efficiency detectors and dark counts, the density matrices for the other detection results were 
multiplied by constant factors, and added to the density matrix for the desired detection result. It was assumed that 
the inefficient detectors register single-photon, two-photon and three-photon states as vacuum with probabilities of 
10%, 1% and 0.1%, respectively. For the detector on mode 2 it was assumed that a single-photon state is registered as 
two or more photons with 0.1% probability, and the vacuum state is registered as two or more photons with 0.0001% 
probability (due to the lower probability of simultaneous dark counts). The appropriate equations to use for the case 
with multiphoton inputs are derived in Sec. 

For perfect sources and detectors, the final probability for a single photon is above the initial probability of 20% 
when the probability for obtaining the detection results is below about 0.7%. Thus, in order to obtain the desired 
detection results, the experiment needs to be repeated roughly 200 times, which is not unreasonable. If we consider 
a final detector that can not distinguish between two, three or four photons, the results are almost identical, so this 
problem is relatively trivial. 

Even photodetectors with finite efficiency do not greatly affect the results. If the first two photodetectors have 
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FIG. 4: The final probability for a single photon versus the probability for obtaining the appropriate detection results for a 
four mode interferometer and p max = 0.2. The black line is the ideal case, the light blue line is that for the case where detection 
results of D = 3 and D = 4 are also allowed, the green line is that if, in addition, the photodetectors have 90% efficiency. The 
dark blue line is that including all of these experimental limitations, plus 0.1% chance of dark counts at detector 2, and the 
red line is that taking into account these experimental limitations, as well as allowing 0.1% probability for two photons in the 
inputs (without dark counts). 

90% efficiency, the final probability for a single photon is reduced by about 0.3%. The greatest problems are dark 
counts at detector 2, and multiphoton components at the inputs. If a dark count rate of 0.1% is allowed, then the 
maximum single-photon probability is reduced below 0.23%. For small values of e the single-photon probability drops 
dramatically, rather than approaching the maximum value. The results are similar if a two-photon probability of 
0.1% is allowed in the inputs (while the vacuum probability is decreased by 0.1%). The single-photon probability 
again drops for small values of e, and the maximum single-photon probability is less than 0.22%. For two-photon 
probabilities of 0.4% or more it is not possible to obtain any increase in the single-photon probability . Thus we see 
that the main problems with experimental realisations will be two-photon components in the input and dark counts 
at detector 2. 



VIII. NO-GO THEOREMS 

In this section, we prove a number of no-go theorems for post-processing via linear optics and photodetection. Note 
that one limitation of the scheme given in Sec. I VII is that it only gives improvements in R for four or modes. In fact, 
it is impossible to obtain improvements for fewer than four modes [m]. This result may be shown in the following 
way. First consider the case D — 0. Then there is only one term in the sum for cq, and Co = K'Pq. The expression 
for ci becomes 

N 

N 

<X'i? in ^F |A lfc | 2 

fe=i 
= K'R in P 

= c i?in. (44) 

Thus we have shown that i? ut < Rin, so c\ < p max - Hence there can be no improvement in the photon statistics if 
zero photons are detected. 

This result can also be shown in a more intuitive way as follows. First note that an arbitrary U(7V) interferometer 
can be obtained using a line of N — 1 beam splitters followed by a \J(N — 1) interferometer (Fig. [SJ). This is an 
immediate consequence of the algorithmic construction of arbitrary U(iV) interferometers from beam splitters [1^. If 
the N — 1 modes upon which the U(iV — 1) interferometer acts are those that are measured, then we may omit the 
\J(N — 1) interferometer entirely (because detecting zero photons at the output of this interferometer is identical to 
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FIG. 5: A U(iV) interferometer can be represented by a U(iV — 1) interferometer preceded by N — 1 beam splitters. This figure 
shows the example for N = 5. 



detecting zero photons at the input). Thus this case may be reduced to the case of a line of beam splitters where zero 
photons are detected at each stage. 

The case of a line of beam splitters with vacuum detection may be deduced from the case for a single beam splitter. 
As was shown above, with a single beam splitter there is no improvement in the ratio between the probabilities for 
detecting one and zero photons. It is easily seen that the same result holds if there are nonzero probabilities for 
photon numbers larger than 1 in the inputs (corresponding to photon numbers larger than one in the output). 

Thus, if we have a line of beam splitters, the ratio between the probabilities for one and zero photons in the output 
can not be increased above the maximum of that for the inputs. This result implies that the probability of one photon 
in the output can never exceed p ma x- This result holds for a line of beam splitters, and therefore for an arbitrary 
U(iV) interferometer. 

We can also obtain a similar result for the case D — 1, provided all the input pi arc equal. If the single photon is 
detected in mode to, then 

c =K'J2 i P " laX p o lAmfel 2 = K'R m P . (45) 

K 

The value of c\ is given by 

k l-l^k 

< ^K'Rf n P £ |Ai;A mfc + A lfe A m/ | 2 

k,l 

= \k'rIp ]T [|A lz | 2 |A mfc | 2 + |A lfe | 2 |A m ,| 2 + A lz A mfe A lfc A^ + A^A^ fe A lfe A m/ ] 



k.l 



= \k'rIp 



£ | A 1Z | 2 £ \^k\ 2 + £ | A lfe | 2 ]T | A ml | 2 + ]T Au A; nl A mkA* lk + A U A ml £ A *nk A U 

K'R? n P . (46) 



In the last line we have used the fact that and A m k are orthonormal. Thus we again find i? ou t < R'm, so c\ < p max . 

These results can be used for an alternative proof that no improvement is possible for the case of a single beam 
splitter. We have shown that detecting zero photons does not give an improvement, and if one photon is detected, 
then we must have M — D = 1 or 0, so there again can be no improvement. 

We can also eliminate the case of a three-mode interferometer, though the reasoning is not as straightforward. First 
note that an input with probability pi of a photon can be obtained by randomly selecting between a source with 
efficiency p max and the vacuum. That is, with probability q — Pi/p m&K we use the source with efficiency p maX 7 an d 
with probability 1 — q we use the vacuum state. If we discard the information about which source was used, this 
is obviously equivalent to a source with efficiency pi. Hence the value of c\ for the source with efficiency pi is the 
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weighted average of the values of c\ for the cases where the efficiency is p max and zero. Thus the maximum c\ must 
be obtained with all of the nonzero Pi equal to p max . 

Therefore, in considering the three-mode interferometer, we can let all the nonzero pi be equal to p max - If ah the 
Pi are nonzero, then we can use the result showing that there is no improvement with one photon detected and all pi 
equal. The other alternatives for detection have D — or M — D < 1, so there can be no improvements in these cases 
either. If one or more of the pi are zero, then the only detection alternatives are with D = or M — D < 1. Thus we 
have shown that there can never be an increase in the probability of a single photon if less than four modes are used. 



IX. UNSOLVED PROBLEMS 



Although it seems that we were able to answer most of the relevant questions concerning the possibility of 
improving the efficiency of single-photon sources, there are, in fact, still a number of open questions we would now 
like to address. The two main unsolved problems for this post-processing are: 

1. Is it possible to increase the probability for a single photon, regardless of the value of p max ? 

2. Can the single- photon probability be increased without adding a multiphoton component? 

At this time the indications are that the answer to both these questions is no. We have performed numerical 
searches for interferometers that give improvements for p > 1/2. These searches have been unsuccessful, indicating 
that it is not possible to obtain an improvement for p > 1/2. We have not been able to prove this assertion; however, 
we can show that there are various implications if there is any value of p max such that it is impossible to obtain an 
improvement. 

First note that it is sufficient to use pi = p max in the input modes. As discussed above, for a given interferometer 
the maximum improvement will always be obtained with all of the nonzero pi equal to p max . It is possible to obtain a 
vacuum state from inputs with efficiency p max , simply by using a beam splitter and conditioning on detection of two 
photons at one of the outputs. Therefore, if there is an interferometer that achieves a certain result using inputs with 
Pi = or p max , there will always be another (expanded) interferometer that achieves the same result with pi = p max . 

Using this simplification, the expression for c ni simplifies to c ni — K"d ni R™^ /nil, where 

d ni = \ S *>™\ 2 - ( 4? ) 

s:S s =S„ 

The values of d ni are independent of p max , and depend only on the interferometer and combination of detection 
results. There is an improvement in the probability of a single photon if 

N 

dxXk+YjdniK/nil. (48) 

rii=2 

Let po be a value of p max such that there is an improvement in the probability of a single photon, and let the 
corresponding value of i?; n be Rq . Then there exists an interferometer and combination of detection results such that 



N 

dx >d + Y, d nx R^/m\. (49) 

rii=2 

Since each of the d n are positive, and the right-hand side is increasing as a function of _R; n , we find that (|48l) is 
satisfied for all < p max < po- Thus we find that, for any value of p max such that there is an improvement, there is an 
improvement for all smaller values of p max . In turn this result implies that, if there is no improvement for p max = po, 
then there can be no improvement for larger values of p max . 

It is also possible to show that, if it were possible to obtain an improvement with no multiphoton contribution, there 
would be no value of p ma ^[2l) for which we can not obtain an improvement. To show this, note that zero multiphoton 
contribution implies that d ni = for ri\ > 2. Therefore, if this improvement is possible for p max = po, then Eq. I|49|) 
becomes simply 

di > d . (50) 

Similarly, the condition to obtain an improvement for any other value of p max is simply d\ > do, which is automatically 
satisfied. In addition, because d ni = for n\ > 2, c ni = for ni > 2, for any p max . 
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Thus, if it is possible to obtain an improvement for some value of p max while maintaining zero multiphoton contri- 
bution, then it will be possible to obtain an improvement for all values of p max . In addition, the relative improvement 
in R is independent of p max . To see this, note that 

-Rout _ c\/ cp _ <h /g-js 
Rin R'm d() 

A further implication is that it would be possible to obtain an output state that is arbitrarily close to the pure 
single-photon state. As the output from the interferometer has no multiphoton contribution, outputs from N of these 
interferometers may be used as the input to another, thus increasing R by a factor of (di/do) 2 . Further iterations 
may be used to increase R by a factor of di/do to any arbitrary power, thus obtaining a final probability for a single 
photon that is arbitrarily close to 1. 

At this stage there is no known scheme that can give an improvement in the probability for a single photon 
while maintaining zero multiphoton contribution. As shown above, if this were possible for any value of p max , 
then it would be possible for all values p max . As increasing the single-photon probability without maintaining zero 
multiphoton component is a less difficult problem, if it were possible to obtain an improvement while maintaining zero 
multiphoton component, for all values p max there would be an enormous range of schemes that give improvements 
without the constraint on the multiphoton component. Such a wide range of schemes would be relatively easy to 
find numerically; the fact that numerical searches have failed to find any scheme that gives an improvement in the 
single-photon probability for p max > 1/2 therefore implies that it is extremely unlikely that there is any scheme that 
gives an improvement while maintaining zero multiphoton contribution. Nevertheless, these numerical results are not 
sufficient to rule out this possibility. 



X. MULTIPHOTON INPUTS 



The majority of this study is based upon inputs from photon sources that have zero multiphoton contribution. 
It is also possible to derive results for inputs with nonzero probabilities for two or more photons, but this case is 
more difficult. The simplest case is for a beam splitter with multiphoton inputs. Let us denote the probability for m 
photons in input mode i by Pi m . Then the input state may be written as 



P 



(AT) 



y^,pikP2i\ki)(M\ 



k.i 



E 

fej 



^i(at)'=(at)'|00)(00|(a 1 )' £ (« 2 ) i . 



The beam splitter transformation gives 



5 (N) 
rtrans 



E 

k,l 



P\kP2 l , 
kill 



(A n a\ + A 21 <4) fe (Ai2<i} + A 22 a 2 )'|00)(00|(A lia i + A; i a 2 ) fe (A 12 a 1 + A* 22 a 2 ) L 



(52) 



(53) 



Expanding in a series and conditioning upon detection of D photons gives 

j D—m A m A k—D+m^l—m 



JN) 
Pout 



- K^p lk p 2l k\l\D\{k + I - D) 



k.l 



min(D,0 

E 

m— m aX (X) — fc,0) 



a u-m a m A K 
iv l 2 yv 2 l 



(D - m)\m\(k - D + m)\(l - m)\ 



= KD\ 



\^22? D 

|A 12 | 2 ^ 



J2P^P2ikm(k + l-D)\\A 12 \ 2k \A 22 

k,l 
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min(D,;) 

E 

m— max(D- k,0) 



(-1)™|A 12 /A 22 | S 



{D - m)\m\[k - D + m)\{l - m)\ 



, (54) 



where if is a normalisation constant. If K is omitted, the trace gives the probability for this detection result. This 
is the expression used to calculate the numerical results in Sec. I VIII 

The transformation for a multimodc interferometer is also straightforward to determine. For this case the input 
state may be written as 



N 



(g)[5> im HH] 



i—l rn 



Ep n(«!) s *i°)(°in^r i 



(55) 
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where P s — Yii(Pi,si) Si ■ From this point on we use the abbreviated notation P' s = P a j \\ i s,!. Note that may take 
any value > 0, rather than simply and 1. Now applying the interferometer transformation &' i— > A T d^ gives 



r trans / . s 



- " 11 \ . 

k 

After detection on modes 2 to N, the final state obtained is 

N 



n e a 4 loxoin e a ^ 



(56) 



/C ) = E^KXml. (57) 

ni=0 

with 



Yi 1 



h ' E ^l^.nl 2 , (58) 



ml 



where S s ^n = per(A[n, s]) and K' is a normalisation constant. This result is very similar to the case for inputs with 
no multiphoton contribution. The only difference is that values of Si larger than 1 are now permitted, and there is 
the additional dividing factor of Yii Si\. We do not use this result in this paper, but it is a useful general form. 

Only one of the no-go theorems still applies for the case where multiphoton inputs are allowed. For the case of a 
beam splitter, the multiphoton contributions in the input give multiphoton contributions in the output. Therefore, 
if zero photons are detected, there can be no improvement in the ratio between the probability for one photon to 
the probability for zero photons. In the multimode case where zero photons are detected, it is possible to decompose 
the interferometer as in Fig. [31 then omit the U(N — 1) interferometer. At each beam splitter in the chain the 
ratio between the probabilities for one and zero photons is not increased, so the final ratio can not be above the 
maximum for the inputs. Nevertheless, the result in this case is not as strong as in the case where the inputs have no 
multiphoton contribution. Proving that the ratio between the probabilities for one and zero photons has not increased 
does not prove that the absolute probability for a single photon has not increased. The problem is that it is possible, 
in principle, for the multiphoton component to be decreased sufficiently that the probability for a single photon is 
increased. 

If the inputs have multiphoton contributions, it is not impossible to obtain a perfect single-photon output state. 
In particular, consider a state po = J2i <&K)(*I) such that qo = 0, qD+i > 0, and qi = for i > D + 1. If this state 
is combined with the vacuum at a beam splitter, and D photons are detected, then the output state will be a pure 
single-photon state. This result also demonstrates that there is no initial value for the single-photon probability that 
can not be improved upon. 

The drawback to these results is that states such as po would be very difficult to produce in the laboratory It is 
likely that there is some measure of the quality of the state that can not be improved upon using linear optics and 
photodctection. However, it is difficult to determine what measure this would be. For example, it is clear that IT can 
be improved, because po can be heavily super-Poissonian. The same considerations rule out other simple possibilities 
such as the entropy. Finding a measure that is non-increasing under linear optics and photodetection is a promising 
direction for future research. 



XI. CONCLUSIONS 



Triggered single-photon sources produce an incoherent mixture of zero and one photons, with much smaller prob- 
abilities for two or more photons. Provided the multiphoton contributions in the inputs may be ignored, we have 
shown that it is possible to significantly increase the probability for a single photon by using post-processing via linear 
optics and photodetection. This method has the drawback that it produces a significant multiphoton component that 
is comparable to that for the Poisson distribution. 

We have shown that there are severe limitations on what post-processing can be performed. In particular, there is 
an upper limit on the increase in the probability for a single photon. This upper limit can not be achieved, but the 
method we have found for increasing the probability for a single photon gives the same scaling with the number of 
modes. It is likely that this method achieves the maximum increase in the probability for a single photon. This result 
is indicated numerically, but has not been proven. 

In addition, it is impossible to obtain an increase in the probability for a single photon using a single beam splitter. 
Alternatively, if zero photons, or one less than the maximum input photon number are detected, it is again impossible 
to obtain an improvement in the probability for a single photon. In the restricted case that all the inputs are identical, 
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it is impossible to obtain an improvement in the probability for a single photon if one photon is detected. These no-go 
theorems are sufficient to prove that at least a four- mode interferometer is required to obtain an improvement. 

Another important no-go theorem is that it is impossible to obtain a perfect single-photon output with imperfect 
inputs. It must be emphasised that this no-go theorem is only for mixed states with no multiphoton components. 
If we relax these constraints, by considering pure superposition states of zero and one photon, then it is possible 
to obtain a pure single-photon output. Alternatively, some multiphoton states can be processed to yield a perfect 
single-photon output. However, it must be emphasised that it is not likely that pure input states, or the appropriate 
multiphoton states, can be produced experimentally. 

There are also a number of important unsolved problems. It is currently unknown whether there is an upper limit 
(less than 1) to the initial probability for a single photon such that it is possible to obtain an improvement. It is 
also unknown if it is possible to obtain an improvement in the probability for a single photon while maintaining zero 
multiphoton contribution. If such a scheme were possible it would be very significant, because it would be possible 
to obtain a state arbitrarily close to a single photon state from arbitrarily poor input states. 
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